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BALLISTIC  RESEARCH  LABORATORIES 

N-‘ 

'  I 

MEMORABBUM  REPORT  NO.  I506 


ADStiegler/sec/eLh 
Aberdeen  Provl3ag  Ground,  Md. 
Novaaber  I96O 


A  MATHEMATICAL  PQRMULATIOK  FOR  ORDVAC  COMPUTATION 
OF  THE  SINGLE-SHOT  KELL  PROBABILITIES 
OF  A  GENERAL  MISSILE  VERSUS  A  GENERAL  AIRCRAFT 

ABSTRACT 

This  report  is  an  up-dating  of  BRL'  Report  No.  867.  That  report  gives 
a  description,  of  the  physical  model  and  the  random  number  generators  used 
in  this  report. 

The  methods  of  BRL  Report  No.  86?  have  become  inadequate  for  the  so¬ 
lution  of  the  problem  with  all  the  recent  modifications  and  refinements 
and  have  been  superseded  by  the  methods  discussed  in  this  report. 

The  main  body  of  the  report  is  devoted  to  a  mathematical  discussion 
of  the  problem  and  appendices  are  included  for  the  programming.  The 
equations  were  formulated  by  Messrs.  De  Gray  and  Stubbs.  The  following 
assumptions  were  made:  the  missile  and  target  are  traveling  in  straight 
lines,  the  rear  ziape  of  the  cone  fuze  pattern  has  been  caultted,  the 
azimuth  and  elevation  angles  are  not  necessarily  standard  for  all  missiles. 

The  discussion  has  been  limited  to  the  fixed-angle  (cone-type)  guided 
missile  fuze,  but  another  fuze  (the  proximity  fuze)  is  being  analyzed  and 
will  be  reported  upon  in  due  course. 
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INTRODUCTION 


The  problem  is  to  determine  the  probability  of  MU  by  blast  and 
fragmentation  of  a  general  missile  against  a  general  target  with  a 
computing  machine  code  simulating  the  real  situation.  The  blast  Is 
determined  against  the  aircraft  as  a  whole,  while  the  fragmentation  Is 
determined  against  various  vulnerable  components  of  the  aircraft  (e.g., 
the  pilot,  bomb,  etc.).  In  order  to  make  the  computing  machine  code  as 
general  as  possible  four  component  types  have  been  allowed  for  and  space 
has  been  left  for  six  more. 

Essentially,  the  code  is  composed  of  three  general  sections.  The 
first  section  deals  with  reading  and  scaling  the  input  data  and  computing 
the  elevation  and  azimuth  angles  of  the  missile,  with  respect  to  the 
target’s  co-ordinate  system.  The  second  section  computes  the  burst  point 
(l.e.,  the  point  at  which  the  warhead  will  detonate)  and  the  effects  of 
blast.  The  third  section  computes,  the  fragmentation  Mils  on  the  indi- 
vidml  components. 

Besides  the  above  there  are  three  other  sections  which  vary  for  each 
target.  These  are  the  comblnatorials  section,  the  vulnerable  areas  section, 
and  the  "variable  numerators"  section.  Each  of  these  will  be  treated  in 
more  detail  below. 

The  final  section  of  the  code  is  the  output  section,  which  prints  out 
probabilities  of  kill  averaged  for  the  given  number  of  rounds  desired. 

NOTE;  A  pictorial  representation  of  the  mathematical  model  can  be  found 
in  BRl  Report  No.  86T. 

INPUT 

The  input  section  reads  in  the  input  arid  scales  it,  initializes  the 
code,  and  computes  the  elevation  and  azimuth  angles  of  the  missile.  The 
input  consists  of  the  following: 

Shields  (Co-ordinates  of  ellipsoid  center  followed  by  axis  lengths) 
Vulnerable  Component  centers  (Co-ordinates  of  each  center) 
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Glitter  Points  (Co-ordinates  of  each  point) 

Blast  (Ellipsoids,  as  above  for  Shields) 

Number  of  Vulnerable  Components 
Number  of  Shield  Ellipsoids 
Number  of  Blast  Ellipsoids 

Number  of  each  specific  vulnerable  component  (Number  of  Pilots, 
number  of  Bombs,  etc.) 

Fragment  spray  angles  measured  from  the  missile  nose 
Number  of  fragments  per  steradian  in  each  spray 
Initial  fragment  velocity  in  each  spray 
Fuzing  and  guidance  errors  (Standard  Deviations) 

Missile  Speed  (ft/sec.) 

Target  Speed  (ft/sec.) 

Total  number  of  rounds  to  be  computed 

Fuze  half  angle  (only  necessary  for  fixed  angle  fuze) 

Iteration  Center  Co-ordinates 

Ballistic  co-efficients  of  the  fragments  in  each  spreiy 
Fragment  mass  in  each  spray 

Distribution  of  Approaches  (See  Elevation  and  Azimuth  Calculation 
for  full  explanation) 

Homing  Points  (Co-ordinates) 

Delay  Time  (sec.) 

Variable  Numerator  Sentinel  (See  pp  16-I9  for  explanation).  For 
scaling  and  location  of  the  above  quantities  see  Appendix  A.  N.B.  The 
shields,  component  centers,  and  glitter  points  remain  static  for  any  one 
target  and  are  therefore  read  in  only  once.  Also,  all  points  are  given 
vith  respect  to  the  target. 

BIEVATION  AND  AZIMUTH  CAICUUTION 

f 

The  target  is  assumed  to  be  moving  in  a  straight  line  on  the  Y  axis 
In  the  -Y  direction.  The  missile  is  assumed  to  be  approaching  the  target 
on  a  random  straight  line.  The  Elevation  is  defined  as  the  angle  between 
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the  missile  velocity  vector  Vj^  and  the  XY  plane  of  the  target  co¬ 
ordinate  system.  The  center  of  the  co-ordinate  system  is  the  center 
of  gravity  of  the  target  and  the  plus  X  direction  is  to  port  and  the 
plus  Z  direction  aloft.  _ 


An  elevation  angle  between  +  degrees  and  -  degrees  la  computed. 
This  angle  is  computed  as  follows: 

Choose  a  uniform  random  number  0  <  p  <  1,  and  then  form  the  function 
l-2p^,  which’  satisfies 

-1  <  l-2p^  <  +  1  . 

(since  all  angles  are  in  degrees,  the  degree  sign  has  been  omitted  in  the 
following  discussion.)  Multiplying  this  function  by  gives  an  elevation 
angle,  E,  such  that: 

-  E^  <  E  <  + 

The  sine  and  cosine  of  E  are  computed. 

The  Azimuth  angle  A  is  the  angle  between  the  missile  velocity  (Vj^) 
and  the  YZ  plane  of  the  target.  In  the  simulation  these  A’s  are  randomly 
computed. 

Having  a  uniform  random  number  p  such  that 

0  <  p  <  1 

form  Tlf  such  that  \|f  =  l-2p  then  -  1  <  ij?  <  1  uniformly. 

Now,  for  the  first  given  percentage  form 

fj^(p)  =  Aj^'i/ 
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whose  limits  are 

-  A.  <  f.  fo^  <  +A- 

±  ~  X"  '  ~  ± 

For  the  second  percentage  define 

^  <  Oj  (A^  -  A^)  Hr  - 

f2(py  = 

t  >  0-  (A^  -  A^)  lif  +  A^ 

whose  limits  are 

1. 

-  A^  <  fgCp)  ^ 

and 

+  k^<  fgCp)  <  +  Ag 

and  in  general  for  the  ith  percentage  define  fj^(p)  as 

t  <0:  (A^  -  k^_^)  t  - 

fl(p)  = 

i  >  0:  (A^  -  A^_^)  t  +  A^_^ 

provisions  will  he  made  for  seven  (7)  distributions  of  A.  The  A^'s 
and  the  percentages  are  given.  For  1  =  1,  A^  =  0.  Then,  F^(p)  =  A, 
where  A  is  the  azimuth  angle. 

Compute  the  sine  and  cosine  of  A  and  then  go  to  the  fuzing  and  blast 
section. 

NOTE;  The  uniform  random  number  generator  Is  described  in  BEL  Report  No. 

855. 

FUZING  AND  BLAST 

The  following  description  involves  two  co-ordinate  systems,  that  of 
the  target  and  that  of  the  missile.  Tie  first  (X,  Y,  Z)  system  has  its 
origin  at  the  center  of  gravity  of  the  target.  The  second  (Q!^,0I2/0{^) 
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system  has  its  origin  at  the  center  of  gravity  of  the  missile  warhead. 

All  data  points  are  given  in  the  system  (i.e,,  relative  to  the 

target) .  Both  missile  and  aircraft  are  traveling  a  straight  course  with-’ 
out  pitch,  yaw,  or  roll. 

Since  the  solutions  of  the  fuzing  equations  are  more  easily  obtained 
with  respect  to  the  missile,  and  all  points  are  given  with  respect  to  the 
target,  the  points  must  be  mapped  Into  the  missile  co-ordinate  system. 
Direction  cosines  of  the  target,  missile,  and  relative  velocity  vector 
are  required.  The  following  conventions  have  been  adopted  in  symbolizing 
these  direction  cosines: 

The  cosines  of  the  angle  between  the  vector  and 
the  first  axis  of  the  system, 
m:  The  cosine  of  the  angle  between  the  vector  and 
the  second  axis  of  the  system, 
n:  The  cosine  of  the  angle  between  the  vector  and 
the  third  axis  of  the  system. 

Unprlmed:  Direction  cosines  of  the  target. 

Primed:  Direction  cosines  of  the  missile. 

Double  Primed:  Direction  cosines  of  the  relative  velocity  vector. 

A  subscript  defines  whether  any  item  is  with  respect  to  the  target  (T)  or. 
the  missile  (m). 

Thus  there  are  six  sets  of  direction  cosines.  These  are; 

(V  V  "^T^  ^^T'  “t"  ^"^T'  “t" 

(VvV  "m' "S' 

Each  is  fully  defined  by  the  above  conventions. 

Since  the  target  is  assumed  to  be  traveling  in  a  -Y  direction  and 
the  missile  is  assumed  to  be  traveling  in  a  -  ot,  direction, 
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0 


> 


■1 


0 


I 


T 


,  ^  = 


Hi 


=  1 , 


“m 


-  0  , 


Hi 


=  0 


The  direction  cosines  of  the  a^es  with  respect  to  the 

X,Y,Z  axes  are; 


II 

cos  E  sin  A 

=  cos  A 

7.1  = 

=  -  sin  A 

sin  E 

II 

CJ 

-cos  E  cos  A 

pg  =  sin  A 

’'a  = 

=  cos  A  sin  E 

(1) 

II 

sin  E 

T) 

II 

o 

^5  = 

=  cos  E 

Thus,  the 

other  direction  cosines  are: 

-  cos  E  sin  A 

■^1 

=  cos  E  cos  A 

= 

cos  E  cos  A 

“m 

=  - 

sin  A 

(2) 

-  sin  E 

“m 

^  — 

sin  E  cos 

A 

if.  t 

1  >  1 

i> 

j  cos  A  cos 

e" 

cos  E  sin  A 

1  Jr  1 

■^T 

1  >  1 

I^r| 

I^'kI 

1  >  1 

1  >  « 

s;  = 

1  ^tI  ^ 

•w." 

“t 

_  I'^mI 

cos  E  COB  A  + 

hil 

I^hI 

I^rI 

(3) 

1  >  1 

t  ^  J 

^  = 

cos  A  sin  E 

T\” 

_ 

sin  E 

> 

^R 

where 


1  ?T  1  f 

|2 

> 

> 

r 

^rI  V 

|v^l  + 

^mI  -  = 

cos  A  cos  E 

Having  determined  the  above  direction  cosines^  choose  the  first 
"glitter  point" .  A  glitter  point  is  a  point  on  the  target  which  is  a 
good  reflector  of  waves  emitted  by  the  fuze  in  the  missile.  Now  since 
the  glitter  point  is  given  with  respect  to  the  center  of  gravity  of  the 
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target  and  the  missile  actually  ’'homes”  or  alms  at ' a  different  point  (viz, 
the  geometric  center  of  all  the  glitter  points),  subtract  the  coordinates 
of  this  "homing  point"  H  from  those  of  the  glitter  point  GP,  thus  trans-  • 
lating  the  center  of  the  system  from  the  center  of  gravity  of  the  aircraft 
to  this  homing  point.  Thus; 


GP 

^GP 

^GP 

^GP 

(4) 

^GP 

N.B.  This  homing  point  varies  for  each  of  the  three  distributions 
of  A  and  is  selected  at  the  time  at  which  A  is  computed. 

We  next  rotate  this  glitter  point  into  the  missile  co-ordinate  system 

; 

I  with  its  origin  at  the  target  homing  point  by  the  equations: 


^P  + 

^2  ^GP 

^5  ^GP 

^GP 

Pi 

^GP  ^ 

^2  ^GP 

^5  ^GP 

^GP 

’'l 

^2  ^GP 

^5  ^GP 

"  ^GP 

where  the  X’s,  p's,  and  7's  are  defined  by  equations  (l)  above. 

The  path  of  H  (the  homing  point)  goes  through  a  point  (o,  o^,  Q(^)  when 
it  crosses  the  plane.  The  path  of  GP  (the  glitter  point)  goes  ttoough 

a  corresponding  (^Qpi  ^  ^pGp'  ^  second  point  is  not 

necessarily  on  the  plane.  The  path  of  H  does  not  pass  throu^  (0,0,0) 

because  of  the  guidance  errors  and  which  are  chosen  randomly  from  a 
normal  distribution.  Translating  to  a  co-ordinate  system  whose  center  Is  at 
the  missile  warhead  center  of  gravity  gives 

(a)  ^  ‘^GP 

(b)  Ogjjp  “  ^  ^GP 

(c)  OijQp  =5:3+  ZiOtjjjp 
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where  missile  warhead  ceatfeercof  ^igiiavlty  In  the 

homing  point  syst^. 

The  equation  of  the  line  traversed  by  the  glitter  point  is 

~  ^IGP  _  ?2  “  *^GP  ^  ~  °SgP  /„x 

“iGP  '  ^IGP  ^GP  "  ^^2  SgP  "  ^“3  ^  ^GP^ 

Substituting  equations  (6)  into  equation  (7)  gives: 

'  °1GP  ^  °2  ~  ^GP  ^  '  °^GP 

This  line  has  direction  numbers 
P  = 


Q  = 

R  =  Og  " 

and  direction  cosines  m|^,  defined  by  (5).  Then 

5  =  -!L_  .  _5_ 

A  “2  -  °2  ^  ■  “3 

Solving  for  Og  and  in  terms  of  gives 

m" 

^  =  a«  +  0^^ 

“5  =  “5  ^ 

Hence,  by  virtue  of  6(b)  and  (c): 


m” 

*  '^“2GP 

“3GP 

=  a* 

^  ir 

“1 

* 

(7*2) 


(7-5) 


(8) 
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Substituting  for  in  (8)  by  use  of  (6)  a  gives 

.n*' 

“m 


°i2GP 


“5GP 


=  A  + 


»  B  + 


^6P  ^ 

“iGP  IT, 


(9) 


where 


A  =  o£^  +  ^, 


‘2GP 


■^OP 


B  -  aj  +  *^GP  ■  ^GP 


these  equations  (9)  give  e3q)res8ions  for  the  glitter  point  in  the  a  - 
system  with  origin  at  the  missile  warhead  center  of  gravity. 

Now  in  the  solution  of  the  fixed  angle  guided  missile  fuze,  we  have 
the  problem  of  finding  the  intersection  of  a  line  (the  line  traversed  by 
the  glitter  point)  and  a  cone  (the  fuze  pattern  of  our  fuze).  The  equations 
of  these  are: 


X  -  X’ 


GP 


and 


2^2 

“2  S 


y  -  Y’ 
— — 

“m 


GP 


GP 


2  2 
tan.  Q 


(10) 


(11> 


respectively,  where  0  is  the  fuze  half -angle,  which  is  given.  Here 

0^,  Oj,  are  the  coordinates  of  a  typical  point  on  the  cone  in  a  system 
with  origin  at  the  missile  warhead  c  of  g. 


Since  equations  (8)  are  equations  of  the  line  in  the  missile  co¬ 
ordinate  system  and  the  problem  is  to  find  the  which  satisfies  equation 
(11),  we  merely  substitute  equations  (8)  in  equation  (11 )  and  solve  for 
o^p,  as  follows: 
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^  ^  ^  B  +  ^  ^  (12) 

2  4.  2  ^ 

=  O^Gptan  0 

Sauarlng  and  collecting  ‘tenns  emd  algebraically  manipula'tiTig  giyeg 
an  equation  of  the  form 

ao^Qp  +  2b  Q^jjp  +  c  =  0  ,  (13) 


where 


©lerefore. 


— 5 M — 

m;;  .  nj  -  i;;  tan  0 


°1GP 


2b 


4  ac 


2a 


With  a  =  1  from  equations  (14)  this  becomes 


(14) 


(15) 


OiGp  =  -  b  iyb^  -  c  (16) 

Then  substituting  in  equation  (6) 

“^If  ~  *^0?  "  '^GP 

where  ot^p  Is  the  co-ordinate  of  the  honing  point  when  the  missile 
would  fuze  if  It  only  "saw"  this  glitter  point. 
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We  compute  equations  (4)  through  (17)^  inclusive,  for  all  glitter 
points.  Then,  the  algebraically  smallest  is  chosen  emd  designated 

This  is  the  co-ordinate  of  the  homing  point  at  the  time  when 
the  missile  actually  fuzes.  Then 

:  “is  =  “0.  *  ^  '"e  ®t 


and 

mV 

^  “5  °1B  T! 


(19) 


where  is  a  given  delay  time  in  seconds  and  e  Is  a  nonaally  distributed 
fuzing  error. 


This  gives  the  co-ordinates  of  the  homing  point  at  the  time  of  burst 
in  the  missile  system  with  center  at  the  missile  warhead.  Then 

X*  =  V 

^B  "  /*"  ^  ^2  ^'2  ‘^B^ 

Zj^  =  +  (X^  +  Pg  CKjg)  ^ 

give  the  co-ordinates  of  the  honing  point  in  the  target  co-ordinate  system 
with  center  at  the  missile  warhead.  It  la  necessary  to  have  the  co-ordinates 
of  the  missile  in  the  target  system  with  center  at  the  homing  point.  These 
co-ordinates  are  the  negatives  of  (20). 

But  these  have  the  homing  point  H  as  center  of  the  co-ordinate  system 
and,  therefore,  they  must  be  translated  back  into  the  system  centered  at 
the  center  of  gravity  of  the  aircraft  in  a  manner  effectively  the  reverse  of 
(4).  Hence,  the  burst  co-ordinates  in  this  system  arei 
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h  ’  -  ^  *  h 

N.B»  Equations  (U)  through  (l6)  are  only  pertinent  to  the  fixed  angle 
fuze.  Other  fuzes  would  have  different  solution;  however,  the  remainder 
of  the  equations  (l)  throu^  (21) , are  the  same  for  any  fuze. 

The  effects  of  blast  on  the  aircraft  in  question  can  now  be  considered. 

Before  determining  blast  we  must  first  determine  if  there  are  any 
"variable  numerators".  This  is  done  by  testing  the  "variable  numerator 
sentinel."  If  the  sentinel  is  positive  or  zero,  there  are  no  variable 
numerators,  and  we  determine  blast  effects.  If  the  sentinel  is  negative, 
there  are  variable  nvnnerators,  and  we  find  what  they  are  and  then  determine 
blast  effects. 

Variable  numerators  occvir  in  the  equations  for  the  shielding  and/or 
blast  ellipsoids. 

The  ellipsoids  which  are  tested  in  the  blast  section  define  the  shape 
of  the  aircraft  \mder  study.  The  shielding  ellipsoids  define  the  aircraft 
Itself  while  the  blast  ellipsoids  define  a  volume  such  that  if  the  warhead 
of  the  missile  under  evalmtion  detonates  within  this  volume  the  aircraft 
Will  be  killed.  If  the  warhead  detonates  outside  this  volume  the  aircraft 
will  not  be  killed.  This  is  a  zero-one  probability. 

The  presence  of  one  or,  more  of  the  variables,  X,  T,  or  Z,  in  the 
numerator  of  an  ellipsoid  rotates  the  ellipsoid  so  that  its  axes  are  not 
parallel  to  the  co-ordinate  axes.  It  also  changes  the  lengths  of  the 
semi-axes  from  a,  b  and  c. 

Having  done  the  above  we  go  on  to  determine  whether  there  is  any 
external  blast.  This  is  done  by , substituting  the  burst  point  (Xg,Yg,Zg) 
in  the  equations  of  each  blast  ellipsoid  until 
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for  some  ellipsoid  or  all  ellipsoids  have  been  tested.  If  there  Is 
external  blast,  we  Investigate  internal  blast.  Hils  is  done  in  the  same 
manner  except  that  we  use  the  shielding  ellipsoids  in  equation  (22). 

If  there  is  no  external  blast  or  no  internal  blast  we  test  for 
perforation  or  debris  type  kills.  A  perforation  type  kill  is  defined  as 
occurring  if  the  missile  enters  and  leaves  some  part  of  the  target  and 
explodes  after  exiting,  i.e.,  pictorially. 


Shielding  Ellipsoid 


A  debris  type  kill  is  defined  as  a  kill  resulting  if  and  when  the 
missile  trajectory  intersects  some  part  of  the  target  and  the  explosion 
occurs  before  entry.  Pictorially, 


After  fully  determining  the  above  we  go  on  to  the  fragmentation 
section.  Tills  section  is  Independent  of  the  fuzing  section  and  the  only 
requirements  are  that  we  have  a  burst  point  (Xg,  Y^,  Zg)  in  specified 
locations  and  we  have  the  direction  cosines  given  in  equations  (l),  (2), 
and  (3). 

NOTE;  A  full  explanation  of  the  Perforation  and  Debris  calculations  will 
be  found  in  Appendix  B.  The  normal  random  number  generator  is  described 
in  BEL  Report  867- 
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FMOMENTA-TION  SECTION 


Since  target  travel  must  "be  taken  Into  account^  after  determining 
the  hurst  point  we  incorporate  a  bias  dependent  upon  fragment  speed.  ' 

An  iteration  center  (usually  the  center  of  gravity  of  the  target  (0,0,0) 
but  not  necessarily)  is  given,  and  the  angle  p  between  the  missile  tra¬ 
jectory  and  the  path  the  fragment  would  have  to  travel  if  it  were  to 
strike  this  Iteration  center 


is  determined  from 


-  S'>‘'  *  K  -  h'!'" 


Since 

(x„  - *  (Ip  - \f  *  (Zr  - 

<^E  -  V  _  ■ 

(x^  - *  (Ir  - 

and 

-  \\  _ 

(Xj  -  (2,  -  Zg)" 


are  the  direction  cosines 


of  line  between  (Xp,Yj^,Zj^) . 


Since  we  are  going  to  iterate  Y^,  we  form  the  following; 

B 

S  ■  '’S,  -  ^  -  V® 


(25) 


(24) 


The  Iteration  is  done  only  in  the  Y  direction  since  the  target  is  moving 
only  in  that  direction.  Then  let(Y^  =  Yt)  and  fom 

13  13 
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Cos  p 


(25) 


By  testing  cos  p  against  the  spray  angles,  we  detenulne  which  spray 
the  fragment  would  be  in  if  It^  struck  the  target  at  this  angle  and  we 
choose  the  corresponding  a  (ballistic  coefficient)  and  (initial 
fragment  velocity  at  the  missile  warhead) . 

The  time  it  would  take  the  fragment  .to  traverse  this  path  is  given 
by; 

+ 

t  .  ^  y  ,  (26) 

o 

where 

"v  '  V' -  X^f  +  (Yj  -  1|)^  t  (i^  - 

The  distance  traveled  by  the  target  in  this  time  is  where 

is  the  target  velocity  which  is  given.  The  new  co-ordinate  is  then 

Y*  ^  (2?) 

Where  y_  is  the  original  Y^  of  the  burst  point  computed  in  the  fuzing 
section,  ' 

■  '  We  perform  tMs  iteration  tliree  times  and  now  have  a  new  Y  co-ordl2iate 
subsequent  calculations  in  the  fragmentation  section. 


of  our  burst  point.  This  new  buret  point  (Xg, 


is  then  used  in  all 


We  are  now  prepared  to  treat  each  of  the  vulnerable  components  of  the' 
aircraft.  We  first  find  the  cosine  of  the  angle  between  the  missile  path 
and  the  path  traveled  by  the  fragment  to  the  first  vulnerable  component 
from  an. equation  much  like  (25)  above  except  that  now  we  replace  the 
iteration  center  with  the  vulnerable . component  center  (X^,  Y^,  Z^)  and  we 
have 
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CCS  (Vj„V^) 


V(Xj  -  X^f  *  (Yj  -  Y|)^  (Z^  -  Zgf 


m 


The  denominator  of  this  fraction  is  designated  as  D  and  used  in  later 
calculations.  Again  we  determine  into  which  spray  this  angle  will  fall 
by  testing  it  against  the  cosines  of  the  given  spray  angles,  y^,  y^,. 

We  choose  the  a  and  and  this  time  N  and  m  corresponding  to  the 
spray  into  which  the  angle  falls.  N  is  the  number  of  fragments  per  steradian 
and  m  is  the  square  root  of  the  mass  of  each  fragment.  At  present  the  code 
distinguishes  five  sprays,  thusly: 


Each  of  these  zones  (l,  II,  III,  TV,  V)  has  a  set  of  numbers  associ¬ 
ated  with  it,  that  is,  an  N,  a,  Vq,  and  m.  The  only  restrictions  on  the 
7j^'s  are  that  y^  >  90°,  y^  S  90°  and  any  7^^^  >  7^. 

We  then  check  to  see  if  the  component  is  shielded.  We  do  this  by 
determining  if  one  of  the  shield  ellipsoids  falls  between  the  two  points 
(Xg,  y*,  Zg)  and  (X^,  Y^,  Z^),  where  the  latter  is  the  vulnerable  component 
center  and  the  former  is  the  burst  point.  The  logic  of  this  routing  is 
explained  in  Appendix  B.  These  are  the  same  as  those  used  in  equation  (28). 
If  a  shield  falls  between  these  two  points,  we  say  the  component  is  shielded, 
store  a  zero  in  the  probability  of  killing  that  peirticular  component,  and 
then  go  on  to  studying  the  next  component.  If,  however,  no  shield  falls 
between  these  two  points  the  component  is  unshielded  and  we  determine  the 
probability  of  killing  it  in  the,  following  manner. 
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We  first  compute  the  striking  velocity  of  the  fragment.  This  is 
done  by  computing  the  remaining  velocity  of  the  fragment  at  the  component 
(V),  it  is  given  by  the  formula  ■ 

V  =  Vge-®  .  '  (29) 

Then  the  striking  velocity  is 

-2  2  -  ('‘t  ' 

Vs  °  ''  D  ■.  (30) 

In  order  to  find  the  vulnerable  area  of  the  component  we  must  form  w, 
where  ' 

w  =  Vgm^/^  X  10"*^  (31) 

We  compute  a  given  function  for  the  particxilar  vulnerable  component 
under  consideration,  viz., 

=  f^(v)  '  ■  (52) 

where  is  the  vulnerable  area  of  the  component.  ■ 

Equation  (32)  varies  for  each  vulnerable  component  and  each  target. 

The  equations  describing  the  A^' s  are  given  for  each  target.  .  We  are  now  in 
a  position  to  compute  the  probability  of  killing  this  component .  We  call 
this  quantity  P  and  say  that 

D  is  defined  by  equation  (28),  A  by  equation  (32),  and  N  is  given.  If 
NA  ^  ■ 

the  quotient  — ^  is  greater  than  l6  we  store  a  one  in  the  probability  of  kill 
D 

for  this  component  since  the  function  of  e  here  calculated  would  be  too  ^ 
to  meike  a  significant  difference  in  equation  (33). 


We  then  go  back  and  determine  the  same  things  from  equations  (28) 
throu^  equations  ( 35)  for  all  components  \intil  we  have  found  Pj^'  s  for 
all  vulnerable  components. 

COMBINATORIAIfi 

When  the  blast  kills  on  the  overall  target  structvire  and  fragmentation 
kills  for  all  components  have  been  determined  and  recorded^  we  go  to  the 
section  of  the  code  which  combines  them  Into  various  combinations  of  kills, 
e.g,,  is  there  a  blast  and  pilot  kill? 

We  sum  these  combined  probabilities  for  each  rovmd  and  then  go  to  .* 
the  output  section. 

OUTPUT 

In  this  section  we  determine  if  we  have  computed  the  given  number  of 
roxinds  (B*),  usually  200,  and  if  not  we  go  back  to  finding  a  new  Azimuth 
and  Elevation  and  proceed  throu^  the  entire  code  once  more. 

If,  however,  we  have  computed  kills  for  all  rounds  we  divide  the  summed 
combinatorialB  by  the  number  of  rounds  done  and  print  out  these  average  kill 
probabilities.  We  then  return  to  the  input  section,  where  we  read  a  new  set 
of  blast  elllpsoidB  and  the  other  data  which  follow  these  up  to  and  Including 
the  delay  time  (D^),  and  proceed  as  before  until  finished. 

This  machine  simulation  has  effectively  calculated  the  single  shot  kill 
probabilities  of  various  missiles  against  various  aircraft. 
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APPEaiDIX,A 

-  '  ‘  Is  . 

location  COISTEMTS  SGAT.TNG  MACHINE  SCAUNO 


.  505 

■  ■■■'■  s’  .  / 

•  10"^ 

;•  ^-19 

,  10“^ 

.■  2-19 

505  , 

■B  ■  . 

10"^ 

'  ,  2-^9 

506 

10"^ 

V  ^  2-19' 

.507 

P 

10"^  • 

2-19 

508 

E 

10"^ 

2’^9 

509 

/b  ,  , 

10’^. 

•  :  2^^ 

50K 

F  --  ' 

10“^  : 

0-20 

NOO  -  NTL 

'  Shield  Numerators 

10 

2 

Shield  Denominator 

10”^ 

'  '  2-8 

„-2a 

N80  -  Nil 

Target  Centers 

10 

2  - 

JOO  -  J7L 

Glitter  Points 

lo"^ 

2-29 

J80  -  Jll 

Blast  Numerators 

10"^ 

2-=S 

Blast  Denominators 

10"^ 

2-8. 

FOO  -  FO)+ 

10“^ 

2-15 

FO5  -  FO9 

10'^ 

2-15 

FOK 

a 

X 

10"^ 

2-“ 

0-11 

FOS 

a 

10  ^ 

2 

y 

,,.-5 

0-11 

FQN 

Op 

10 

2 

FOJ 

V 

m 

.  10'^ 

2-« 

Fd" 

10'^ 

2-« 

' 

o-lO 

FOL 

B* 

10 

2 
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APPENDIX  A  (Continued) 


LOCATION 

UUNTiSCri'U 

mruT  suauLNu 

MAiOiliJNJii  DUUiUNU 

FID 

tan^ 

10“^ 

2-fi' 

PU  -  F15 

hi*  h*  ^ 

10"^ 

2-20 

F14  -  PIT 

7i  -  74 

10”^ 

7j^/360  X  2° 

FIB  -  FIN 

“1  "  **3 

10° 

2° 

FU  -  F21 

10'^  .J 

2° 

F22  -  F25 

Percentages 

10° 

2° 

F24 

Variable  Num* 
Bent. 

Positive  - 
No  Variable 

Negative  - 

Variable  Numerators  at 

920 

F25  -  F2J 

Bbmlhg  Points 

10"^ 

2-20 

F30 

D^  In  seconds 

10° 

2° 
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LOCATION 


COHTEHTS 


4N0 

06  008  00  303 

J4NI 

06  000  00  NOO 

ltfI2 

08  000  00  n8o 

4N5 

08  000  00  JOO 

4n4 

08  000  00  j8o 

i|-N5  ' 

08  051  00  poo 

4n6 

03  33J  80  380 

kN7 

1 

00  69L  01  69-1 

4n8 

03  33N  00  36N 

4N9 

00  69L  01  6kw 

4NX 

03  33W  00  360 

4NS 

00  734  00  735 

I4NN 

00  005  00  000 

4NJ  , 

03  3^  00  545 

4NF 

00  6k9  00  6k4 

4nl 

00  F OK  00  300 

4J0 

08  008  00  590 

4J1 

2-22 

4J2 

2-15 

4J5 

10^  X 

4J4 

10^  X  2"^ 

105 

10^  X  2"^^ 

4j6 

10  X  2"^^ 

-  DESCRIPTION 
C.W.  No.  1 
C.V.  N6*  2 
C.V.  no.  ? 

C.w.  Ho,  k 
C.W.  No.  5 
C.W.  No.  6 
BKW.  No.  1 
BKW.  No,  2 
IBKW  No.  1 
IBKW  No.  2 
SKW  No,  1 
BKW  No.  2 
K  for  iteration 
PKW  No.  1 
PKW  No,  2 
C.W.  NRNSR 
KWO 
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LOCATION 


CONTENTS 


DESCRIPTION 


4J7 

560  X  10"^ 

4J8 

12  X  2"^^ 

4J9 

5^5  X  2"59 

UK 

45/560 

US 

90/360 

UN 

180/360 

Uj 

16  X  2"^° 

UF 

Not  Used 

UL 

16  X  2"^^ 

1^0 

00  NOO  00  NOO 

Shield  Numerator  Address 

Ul 

00  000  00  564 

Retvini  Address  Charge 

U2 

00  J80  00  J80 

Blast  Address 

Uj 

00  390  00  390 

Address  of  Z  's 

4Flf 

00  N80  00  n8o 

Target  Center  Address 

00  380  00  380 

P  Address 

0 

4f6 

i 

GP 

Counter  for  glitter  pts. 

4F7 

00  KOO  00  KOO 

Pilot  A  Address 

V 

U8 

00  K40  00  K40 

Engine  A^  Address 

4F9 

00  K80  00  k8o 

Bomh  A  Address 

V 

4FK 

00  KNO  00  KNO 

Fuel  Line  A^  Address 

4PS 

C 

Counter  for  iteration 
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LOCATION 

CONTENTS 

DESCRIPTION  ’ 

300 

“2 

Normal  Random  Number  x  a 

301 

“5 

Normal  Random  Number  x  a 

.y 

302 

6 

i 

Normal.  Random  Number  x  0^ 

303 

s 

Hiamber  of  Shields  x  2~^^ 

304 

T 

-IQ 

BTamiber  of  Targets  x  2  ■ 

505 

B 

Number  of  Blasts  x  2~^^ 

306 

G 

Number  of  Glitter  Pts.  x 

307 

P  . 

Number  of  Pilots  x  2“^^ 

308 

E 

Number  of  Engines  x 

509 

Number  of  Bombs  x  2"^^ 

30K 

F 

Number  of  Fuel  Lines  x  2~^^ 

50. 

3B 

3  times  the  nxunber  of  Blasts 
x  2“^^ 

3ON' 

Po 

Uniform  Random  Number 

3OJ 

nJ 

B*  times  first  Percentage 
x2-^^ 

3aF 

n| 

N?  +  B*  times  second  Percentag 

o-lO 

X  2 

30L 

b 

Present  round  x 

310 

a 

Present  drag  coefficient 

X  2° 

311 

j 

Present  fragment  velocity 

^  x2-^5 

312 

N 

Present  number  of  fragments 

0-15 

X  2 

313 

Present  fragment  weight 

X  10’^  X  2° 

M 


LOCATION 


CONTENTS 


DESCRIPTION 


51*^ 

A 

Azimuth  Angle  divided  by 

560  X  2° 

515 

E 

Elevation  Angle  divided  by 

360  X  2° 

516 

517 

^2 

Direction  cosines  of  Of, 

-1 

axis  X  2 

518 

S 

519 

Pi 

51K 

Pg 

Direction  cosines  of 
'  -1 

axis  X  2 

51 

51N 

’'1 

JU 

’'a 

Direction  cosines  of  a. 

-1  ^ 

axis  X  2 

5IP 

^5 

.1 

31L 

320 

“m 

Direction  cosines  of  Teurget 

with  respect  to  missile 

X  2'^ 

321 

“m 

322 

'ii 

Direction  cosines  of  Relative 

Velocity 

323 

■ft 

Vector  with  respect  to  missile 

X  2"^ 

yzk 

ft' 

as 


;LOCA!EION 

CONTENTS 

DESCRIPTION 

525 

't 

Direction  cosines  of  missile 

526 

“t 

vith  respect  to  target  x  2”^, 

527 

328 

Coordinates  of  present 

329 

^GP 

Coordinates  of  present  glitter 

point  X  2~  (with  respect  to 
target) 

32K 

^GP 

32s 

^IGP 

Coordinates  of  present  glitter 
-20 

.point  X  2”  (with  respect  to 

32N 

'^2GP 

missile) 

32J 

^3GP 

32F 

Y* 

Iterate  burst  point  for 

fragmentation 

32L 

-  coordinate  of  fuze  pt.. 

largest  so  far 

330 

331 

“2B 

Coordinates  of  Burst  Point 

—20 

with  respect  to  missile  x  2“ 

332 

«5B^ 

333 

Coordinates  of  Burst  Point 

—20 

with  respect  to  target  x  2 

33^ 

335 

'\ 

336 

\  -15 

Relative  Velocity  Vector  x  2 

.  \ 
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LOCATION 


CONTENTS 


DESCRIPTION 


557 

A 

(“2  ^2GP  -  ^10P> 

X  2-2° 

■538 

B 

(“3  ^  ^30P  - 

- 

x2-2° 

539 

b 

i''(ni"A  +  n"B)  „-10 

+  n"^  -  tan^  0 

35K 

c 

r^CA^  -H  b2)  „.io 

in''2  t  n’'2  .  Q 

33S 

^2 

+  (h  ■  hH^ 

X  2-^5 

55N 

((Xj  -  +  (h  -  ==  2-^5 

35J 

-  IK  ^2  =  2-^5 

55F 

^1 

-  Y|)^  +  A^  X  2'^^ 

33L 

M 

Number  of  targets  done  so 

far  X  2'^^ 

w 

(V  m  X  10"^)  X  2'^° 

S 

541 

A 

V 

-IS 

Present  A  x  2 

V 

342 

"•  “tW  -  ^B^  + 

^tW  -  *  2“^ 

343 

D 

( (^  -  Xg)  ■^'  -  y^)  + 

X  2'^° 
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LOCATION 


CONTENTS 


DESCRIPTION 


344 

cos  V^  V^ 

3^5 

346 

347 

°T 

350 

TEMPORARY 

353 

STORAGE 

360 

361 

^T 

562 

371 

(Cos  X  2"^. 

372 

(Cos  X  2"^ 

373 

(Cos  7^)  X  2~^ 

374 

(Cos  7^)  X  2'^ 

390  -  593 


580  -  58s 


(Cj^/d)  X  2 

Direction  cosine  of  Relative 
Velocity 

Vector  with  respect  to  the  target 


Present  target  coordinates 

X2-20 
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APPENDIX  B 


PERFORATION  OR, DEBRIS: 

The  problem  Is  to  determine  whether  the  trajectory  of  the  missile 
intersects  one  of  the  target  coinponents  (given  by  the  shielding  ellipsoids) 
and  whether  burst  occurs  before  entry  or  after  exit. 

The  problem  mathematically  is  to  find  the  solution  of  a  line 
intersecting  an  ellipsoid.  The  equations  are 

The  line: 

(a)  X  =  Xg  +  t  ij,  y  =  yg  +  t  m"^,  z  =  +  t  n"^  , 

where  0  t  1  (B,l) 

The  ellipsoid: 


where  and  are  the  direction  cosines  of  the  relative  velocity 
vector  with  respect  to  the  target,  and  the  ellipsoids  are  the  shielding 
ellipsoids  which  are  given. 

Substituting  (a)  (b.I)  into  (b)  (B,1)  gives 


55 


Squaring  the  terms  In  equation  (B.2)  and  combining  terms  gives 


Letting 

A  =  (a^,  ag,  a^)  (B,4) 

®  ®  S2>  Sj) 


where 


equation  (B.?)  becomes 

[a  1  ijt®  +  2[a  .  ojt  +[q  .  o]  -  1  '  a  0  (B.6) 


If  equation  (B.6)  has  two  read,  and  distinct  roots  there  is  either 
perforation  or  debris.  Since  perforation  Is  defined  as  occurring  when 
the  missile  has  entered  and  left  some  part  of  the  target  and  esqploded 


after  exit,  both  roots  must  be  negative,  or  the  following  criteria  ....  ; 
must  hold: 

(a)  (A  •  A)  (G  •  G  -  1)  -  (A  •  G)  <  0 

(b)  ^  •  G  -  1  >  0  (B.7) 

Condition  (a)  (B.7)  meajis  that  there  are  two  real,  vtnequal  roots. 
Condition  (b)  (B.7)  means  that  the  roots  have  the  same  sign  since  A  .A 
must  be  positive.  Condition  (c)  (B.7)  means  that  at  least  one  root  is 
negative . 

Debris  is  defined  as  the  condition  arising  when  the  missile  trajectory 
Intersects  the  target  and  the  missile  explodes  before  entry.  Therefore, 
both  roots  of  equation  (B.6)  must  be  real  and  positive.  Conditions  (a) 
(B»7)  and  (b)  (B.7)  above  must  again  hold.  However,  in  place  of  condition 

(c)  (B.7)  the  following  must  hold 

‘V..  , 

(d)  i  .  G  <  0  (B.7) 

This  means  that  at  least  one  root  is  positive,  since  A  .  A  is 
always  positive. 

If -condition  (a)  (B.7)  does  not  hold  there  is,  obviously,  neither 
type  of  kill. 

This  analysis  is  for  all  ellipsoids  until  conditions  (a)  (B.7) 
and  (b)  (B,7)  and  (c)  (B.J)  or  (d)  (B.7)  are  met  or  sill  ellipsoids 
have  been  tested. 

i 

■SHIELDING; 

Here  the  problem  is,  as  above,  to  solve  for  the  Intersection 
of  a  line  and  an  ellipsoid.  However,  in  this  caae  the  equation  of 
the  line  is  \  (B.8) 

X  =  Xg  +  t(x^  -  x^),  y  =  yg  +  t(y^  -  y^),  z  =  +  t(z^  -  z^) 

where  0  t  S  1  . 
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Substituting  these  values  In  equation  (b)  (B.I),  gives  an  equation 
fl-t miliar  to  equation  (B,6),  The  equation  is 


1  2  r  > 

>  1 

f  >  >  1 

[b  •  B 

.J  t  +  2  [b 

•  GJ  t 

+  |_  G  •  G  J 

(B.9) 


and 

I  =  (b^,  b^,  -b^)  ^  (B.IO) 

and 

(Xj  -  Xg)  (Yj  -  Yj) 

h - - 37-  > 


(B.u) 


The  path  of  the  freigment  must  intersect  the  ellipsoidj  i.e.,  the 
qviadratic  (B.9)  has  real,  roots,  also  the  ellipsoid  must  fall  between 
the  two  points  (l.e.,  neither  point  can  be  within  the  ellipsoid  in 
question).  Therefore,  the  two  roots  t^  and  tg  of  equation  {B.9)  must 
meet  the  following  criteria: 


0  <  ti  tg  <  1 


(B.I2) 


Since  (B,9)  is  a  function  of  t,  and  if  it  is  designated  by  f(t), 
then  the  ellipsoid  shields  the  conponent  if,  and  only  if,  it  satisfies 
the  following  requirements: 


(a) 

(1  .  1)(G  •  ^  -  1)  -  (B  •  G) 

f(0) 

>  0 

(c) 

f(l) 

>  0 

(d) 

f'(0) 

<  0 

(e) 

f'(l) 

>  0 

where  f'(t)  is  the  derivative  of  the  function  with  respect  to  t. 


(B.I5) 


36 


If  condition  (a)  (B.I5)  holds  "both  the  roots  are  real.  If 
conditions  (h)  (B.lj)  and  (c)  (B.13)  hold  then  neither  the  point 
(Xj,  Yg,  Zg)  nor  the  point  (X^,  Zj),  respectively,  are  inside 
the  ellipsoid.  If  conditions  (d)  (B.15)  and  (e)  (B.l 5)  hold  there  is  . 
either  a  jnaxlmum  or  mini  mum  between  0  and  1.  Since  condition  (a)  (B^l^) 
is  tested  first,  this  means  that  there  is  a  root  Between  0  and  1..  Thus 
conditions  (B.lj)  are  synonymous  with  conditions  (B.12), 

■  If  none  of  the  eUlpsolds  meet  the  requirements  (B.15)  the 
component  is  not  shielded.  However,  if  at  least  one  ..ellipsoid 
meets  these  requirements,  then  the  con^tonent  is  shielded. 
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